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ABSTRACT 

We compare high-resolution optical and radio imaging of 12 luminous submillimeter (submm) galaxies 
at a median z = 2.2 ± 0.2 observed with Hubble Space Telescope (HST) and the MERLIN and VLA 
radio interferometers at comparable spatial resolution, ~ 0.3'' (~ 2kpc). The radio emission is used 
as a tracer of the likely far-infrared morphology of these dusty, luminous galaxies. In ~30% of the 
sample the radio emission appears unresolved at this spatial scale, suggesting that the power source is 
compact and may either be an obscured AGN or a compact nuclear starburst. However, in the majority 
of the galaxies, ~70% (8/12), we find that the radio emission is resolved by MERLIN/ VLA on scales 
of ~ l''(~10kpc). For these galaxies we also find that the radio morphologies are often broadly similar 
to their restframe UV emission traced by our i/ST imaging. To assess whether the radio emission may 
be extended on even larger scales, 1'', resolved out by the MERLIN+VLA synthesized images, we 
compare VLA B-array (5'' beam) to VLA A- array (1.5'' beam) fluxes for a sample of 50 /iJy radio 
sources, including 5 submm galaxies. The submm galaxies have comparable fluxes at these resolutions 
and we conclude that the typical radio emitting region in these galaxies are unlikely to be much larger 
than ~ l"(~10kpc). We discuss the probable mechanisms for the extended emission in these galaxies 
and conclude that their luminous radio and submm emission arises from a large, spatially-extended 
starburst. The median star formation rates for these galaxies are ~ 1700 M© yr~^ (M> O.IM©) occuring 
within regions with typical sizes of ~40kpc^, giving a star formation density of 45 M© yr~^ kpc~^. Such 
vigorous and extended starburst appear to be uniquely associated with the submm population. A more 
detailed comparison of the distribution of UV and radio emission in these systems shows that the broad 
similarities on large scales are not carried through to smaller scales, where there is rarely a one-to-one 
correspondance between the structures seen in the two wavebands. We interpret these differences as 
resulting from highly structured internal obscuration within the submm galaxies, suggesting that their 
vigorous activity is producing wind-blown channels through their obscuring dust clouds. If correct this 
underlines the difficulty of using UV morphologies to understand structural properties of this population 
and also may explain the surprising frequency of Lya emission in the spectra of these very dusty galaxies. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation — galaxies: 
starburst 



1. INTRODUCTION 

Since their discovery, luminous submm galaxies (SMGs) 
have been proposed as candidates for the progenitors of 
the most massive spheroids in the local Universe (Small, 
Ivison & Blain 1997; Hughes et al. 1998; Lilly et al. 1999; 
Blain et al. 2002). The recent measurement of the red- 
shift distribution, space densities and clustering of this 
population provides strong support for this proposed rela- 
tionship (Chapman et al. 2003a, 2004; Blain et al. 2004). 
These galaxies have large bolometric luminosities, ~10^^- 
10^^ Lq^ characteristic of ultraluminous infrared galax- 
ies (ULIRGs, Sanders & Mirabel 1996). If their intense 
restframe far-infrared (far-IR) emission arise from dust- 
obscured star formation, then the estimated rates are 
^ 10^ Mq yr~^, sufficient to form the stellar population of 
a massive elliptical galaxy in only a few dynamical times, 
given a sufficient gas reservoir. 

Alternatively, a substantial fraction of the submm emis- 
sion in these galaxies could arise from an obscured AGN 



(e.g. Almaini et al. 1999). It has proved difficult to dis- 
tinguish whether AGN or starburst activity powers the 
dust heating and associated far-IR radiation in these lumi- 
nous submm galaxies (Frayer et al. 1998; Alexander et al. 
2003; Chapman et al. 2003a). Optical and near-infrared 
spectroscopy or X-ray observations have frequently been 
used to search for the signatures of AGN, in both local 
ULIRGs and those at high-redshifts (Sanders & Mirabel 
1996; Fabian et al. 2000; Ivison et al. 2000; Barger et al. 
2001; Frayer et al. 2003; Chapman et al. 2003a; Alexander 
et al. 2003, 2004; Swinbank et al. 2004). However, merely 
identifying the presence of an AGN within a ULIRG does 
not immediately mean that it must be the dominant source 
of far-IR radiation. Energetic arguments must be used 
to estimate what fraction of a ULIRG 's luminosity arises 
from the AGN. 

A much simpler test is available if the far-IR emission 
is resolved - the geometry of the emission from an AGN 
means it is not a natural source to heat dust over an ex- 
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tended region - hence any extended far-IR emission is very 
likely to arise from star formation. Unfortunately, the 
coarse resolution of most far-IR and submm instruments, 
e.g. - 15'' FWHM at 850 /im with SCUBA on the JCMT, 
means that the emission is rarely resolved except in the 
most local galaxies (Le Floc'h et al. 2002). There have 
been recent claims for the detection of submm emission on 
~ 100 kpc scales around some powerful high-redshift AGN 
(Ivison et al. 2000; Stevens et al. 2003, 2004), however, 
these are rare and extreme objects whose characteristics 
may have little bearing on those of typical submm galaxies. 
To disentangle the mechanisms responsible for the far-IR 
emission in the population of submm galaxies at 2: 2-3 
(Chapman et al. 2003a), will likely require sub-arcsecond 
resolution to map emission on kpc scales - well beyond the 
capabilities of current far-IR/submm facilities. 

One way to circumvent the limited spatial resolution of 
far-IR/submm instruments is by exploiting the tight far- 
IR-radio correlation observed for infrared galaxies (e.g. 
Helou et al. 1985; Condon 1992) and the high angular res- 
olution capabilites of long-baseline radio inteferometers, 
such as the Multi-Element Microwave Linked Interferom- 
eter (MERLIN) or the Very Large Array (VLA), to infer 
the sub-arcsecond distribution of far-IR emission within 
submm galaxies. One caveat of this approach is that the 
far-IR-radio correlation has only been demonstrated lo- 
cally on relatively large scales, ~ 50 kpc (Yun, Reddy & 
Condon 2001), and it the precise correlation may break 
down on the smallest scales (M. Yun, in preparation). 
Nevertheless, if a significant extended component of the 
continuum radio emission from the submm population is 
seen then this would provide strong support for the far-IR 
emission being similarly extended. 

The deepest 1.4-GHz observations from the VLA can 
detect the synchrotron emitting disks and nuclear star- 
bursts (formed from the coalescence of radio supernovae 
and their remnants) of a ULIRG such as Arp220 out to 
redshifts of z = 3-4. Indeed, a large fraction of the bright 
submm population at high redshifts are detected as /iJy 
radio sources (Small et al. 2000; Barger, Cowie & Richards 
2000; Chapman et al. 2001; Ivison et al. 2002), as expected 
given their submm luminosities and the local far-IR-radio 
correlation (Chapman et al. 2004). The highest spatial res- 
olution available from the VLA at 1.4 GHz is 1.5'', but by 
combining deep 1.4- GHz observations from MERLIN and 
the VLA it is possible to produce datasets which combine 
both high sensitivity and high spatial resolution, ^ 0.3" 
scales, sufficient to map /iJy radio sources and test the 
extent of the far-IR activity in these galaxies. 

These sub-arcsecond radio maps, which indirectly trace 
the far-IR morphologies, of the SMC population can also 
be compared and contrasted with the restframe UV struc- 
tures visible in HST imaging on similar scales. Such com- 
parisons may help to constrain the extent of obscuration 
in SMGs relative to other samples of high-redshift sources 
selected in the restframe UV (e.g. Adelberger & Steidel 
2000). 

In this paper, we present sensitive MERLIN/VLA radio 



and iJST restframe UV observations of SMGs at compara- 
ble, sub-arcsecond resolution. We discuss the sample and 
observations in §2, describe our main results in §3 and dis- 
cuss our conclusions in §4. We assume a A-CDM cosmol- 
ogy with flo = 0.3, 11a = 0.7 and Hq = 65kms~^ Mpc~^, 
so that larcsec corresponds to 8.4 kpc physical size at 
z = 2.2. 

2. SAMPLE AND OBSERVATIONS 

Our primary observational dataset is the combined 
MERLIN^ and VLA radio map of a 10' x 10^ region cen- 
tered on the Hubble Deep Field (HDF) region, which has 
sufficient sensitivity and resolution to attempt a radio mor- 
phological analysis of the submm galaxies in this area. 

The SMG sample in this region comes from the par- 
ent catalog of submm-detected, optically- faint /iJy radio 
galaxies used in the spectroscopic survey of Chapman et 
al. (2003a, 2004). We identify 14 submm-detected galaxies 
lying within a 7.5'-diameter field centered on 12 36 48.0 
-h62 15 40 (J2000) which have 850/im fluxes brighter than 
4mJy and are detected at 1.4 GHz in the VLA A-array 
observations of this region with a flux of > 40/iJy. This 
radio flux limit should ensure a useful constraint on the 
source morphology from the MERLIN observations. 

2.1. MERLIN Observations 

The deep MERLIN observation of the HDF region 
(Muxlow et al. 1999; Muxlow et al. 2004) comprise ~ 
430 hrs integration at 1.4 GHz of a 10' x 10' field centered 
on the HDF and including the Hubble Flanking Fields 
(HFF). These data were acquired in February 1996 and 
April 1997. The MERLIN data were supplemented with 
42 hrs of 1.4-GHz VLA'' A-array observations (Richards 
2000), and combined and deconvolved in the sky-plane 
due to computational limitations. We use a map with 
a restored 0.3-'' beam which has an rms noise level of 
3.3 /iJy beam" ^. To register the radio and optical data, 
radio sources associated with compact galaxies have been 
used to align the radio map with panoramic ground-based 
imaging (see below). For the HDF, this matching involves 
128 / < 24 optical sources with radio counterparts and 
yields an rms of 0.3'' (Capak et al. 2004). 

2.2. HST imaging 

We next search the HST^ database for deep imaging 
observations of SMGs that lie within MERLIN field. As 
our primary goal is a comparison of the coarse morpholo- 
gies of the sources in the radio and optical wavebands, the 
choice of camera used for the iJST observations is less crit- 
ical than it would be for a detailed morphological analysis 
of SMGs (Chapman et al. 2003b). Hence, we search for 
any observations of SMGs within the HDF MERLIN field 
using the STIS and ACS cameras. 

Due to the intensive study of this field we identify HST 
imaging of 13 SMGs from our sample which lie in the 
MERLIN field and fist these in Table 1. Four of these 
galaxies come from the HST targeted survey of 13 SMGs 

® MERLIN, a UK National Facility operated by the University of Manchester at Jodrell Bank Observatory on behalf of PPARC 

^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 

^The optical data is based around observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science 
Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. 
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Fig. 1. — Comparison of MERLIN/VLA radio data (0.2-0.3^^ synthesized beamsize, shown as contours) for 12 of our SMGs with the HST 
STIS clearSOccd and ACS 6435+2775 imaging (as noted in Table 1). The HST images have been smoothed with a Gaussian of FWHM 0.2" 
to match the effective resolution of the radio map. The first two rows show sources where the radio is clearly extended along a portion of 
the rest-UV imaging, whereas the final four panels on the bottom row reveals point sources in the radio, likely associated with a single UV 
component. Radio contour levels begin at 3a and increase in linear steps scaled to the peak of the radio emission. Each panel is 3^' on a side 
and they are ordered from top-left as in the same sense as Table 1. 



observed with STIS by Chapman et al. (2003b), see Ta- 
ble 1. In addition to these, a further 9 SMGs serendip- 
itously fall in the HDF/HFF region which is covered by 
ACS imaging from the GOODS project (Giavalisco et al. 
2004). 

However, in one of these 13 sources (SMM J 12365 1.76 
+621221.3), the SMC is not uniquely identified, and may 
lie behind a foreground elliptical galaxy (Dunlop et al. 
2004), and an accurate comparison of the morphology us- 
ing the MERLIN/VLA radio and HST in the optical is 
impossible. We exclude this object from further compari- 
son. 

The Cycle 10 STIS imaging of galaxies in our sam- 
ple uses the open filter, clear50ccd (central wavelength 
5 733 A), and have durations of two or three orbits (5.0- 
7.5 ks). The reduction and analysis of these images is de- 
scribed in Chapman et al. (2003b). The resolution of these 
images is 0.06'' FWHM and they have a typical point- 
source sensitivity limit of 27.4 mag. (AB). 

The ACS observations of the SMGs lying within the 
GOODS-N field were obtained from the STScI vl.O release 
(August, 2003), in the ^775 and 6435 bands. The reduction 
of these data is described by Giavalisco et al. (2004) and 
the typical point-source sensitivity is 28 mag. (AB) and 



resolution of 0.07''. 

We register the HST images to the radio coordinate 
frame by aligning them with the deep Suprimecam im- 
ages from Capak et al. (2004) which are tied to the radio 
frame. To achieve this we first smooth the HST images to 
the ground-based seeing, then match all > 5-cr sources (ex- 
cept the SMC), and transform the coordinate grids using 
the IRAF task, geotran. 

The morphological characteristics of this sample have 
been classified by eye from the HST imaging in Fig. 1 
and are presented in Table 1. These rough morphological 
classes are subject to uncertain structured dust extinction 
(Small et al. 1999), and may not represent the true phys- 
ical morphology of the system. 

2.3. Archival observations 

Optical photometry for our SMG sample in the B and R- 
bands (Table 1) was measured from the Subaru Suprime- 
cam imaging published by Capak et al. (2004). We use a 
3" diameter aperture centered on the radio source and the 
limiting magnitudes are B < 26.9 and R < 26.6 (5(j). 

11 of these 12 SMGs in our sample have secure spec- 
troscopic redshifts from the Keck survey of Chapman et 
al. (2004 - Table 1). These enable us to measure physical 
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sizes and luminosities for these galaxies and also calculate 
K-corrections between galaxies in the sample to directly 
compare their restframe properties. UgRIK photometry 
exists for the remaining source (SMM J123646. 1+621449), 
allowing us to estimate a photometric redshift. Using the 
HYPER-z software (Bolzonella et al. 2000), we derive a pho- 
tometric redshift of z = 1.7±0.2. The median redshift of 
the sample is <z>= 2.2zb0.2, representative of that mea- 
sured for a larger spectroscopic samples of submm galaxies 
Chapman et al. (2003a, 2004). 

3. ANALYSIS AND RESULTS 

Using the registered radio and optical imaging we 
show in Fig. 1 the 1.4-GHz MERLIN/VLA contours over- 
layed on the iJ^T images for the 12 SMGs in the joint sam- 
ple. We note that at the median redshift of our sample, 
the HST imaging corresponds to restframe wavelengths of 
1700-1 800 A. The median diameter of the radio emission 
from the 12 galaxies is 0.83 =b 0.14'' or 7.0 ± 1.1 kpc (mea- 
sured above the 3-cr contour. Fig. 1), showing that the typ- 
ical source in our sample is well-resolved at the resolution 
of our MERLIN/VLA map. The radio morphologies for 
these submm galaxies split into two broad classes: those 
dominated by an unresolved component, often centered 
on a sub-component of the i7S' T-optical emission (4/12 
or 33%), and extended structures on scales ~0.5-l''(8/12 
or 67%). Remarkably, the extended radio morphology in 
these latter sources often appears to trace the same UV- 
bright, large-scale structures seen in the HST optical im- 
ages (Fig. 1). 

With sensitive, high-resolution imagery in the radio and 
optical, we can compare the radio emission (as a tracer 
of the dust emission) to restframe-UV emission on kpc- 
scales within the systems. In Fig. 2 we show the radio 
and restframe UV surface brightness profiles along the 
major axes of three representative resolved sources. In 
some cases the UV emission shows broad similarities to 
the radio emission. However, Fig. 2 demonstrates that the 
UV/radio flux ratios can vary significantly over the ex- 
tent of the galaxy, and the SFRs derived from the respec- 
tive wavelengths will differ accordingly. For example, in 
SMM J123621.3+621708 the extended restframe UV emis- 
sion does not trace the more compact radio emission, with 
the peak of the radio emission actually coinciding with a 
clear deficit in the UV emission. Such anti-correlations in 
the radio and UV are similar to those seen in many local 
ULIRGs (Chamandaris et al. 2002). 

To test whether there is any evidence for a cor- 
relation between internal reddening and the distribu- 
tion of obscured star formation traced by the ra- 
dio emission, we construct (6435 — ^775) color maps 
of the 12 SMGs using the ACS imaging of the 
GOODS-N region. Only SMM J123712.0+621325 and 
SMM J123707.2+621408 show redder internal color struc- 
ture which corresponds to the radio morphology. In 
all other cases, the radio emission does not corre- 
spond to any regions with unusually red colors in the 
HST (&435 - ^775) maps. SMM J123707.2+621408 and 
SMM J123712.0+621325 actually show bluer colors in the 
vicinity of the peaks in the radio emission. This suggests 
that the UV emission may not always probe the true site 
of far-IR emission; and indeed the UV-inferred bolomet- 
ric luminosities in SMGs on average underpredict the true 



luminosities by factors ;»10 (Chapman et al. 2004). 
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Fig. 2. — The surface brightness profiles in the radio and op- 
tical wavebands along the major axes of the radio emission in 
three of the resolved SMGs. The dashed line shows the radio 
emission from the MERLIN/VLA map and the solid line shows 
the equivalent distribution of the restframe UV emission from 
the HST imaging. In both cases the surface brightness is aver- 
aged over a 0.3^^ window orthogonal to the major axis. The top 
axis gives the angular size at the mean redshift of these three 
galaxies, the variation in their redshifts produces less than a 
5% change in angular scale between sources. 



A more quantitative test of the degree of obscuration 
within these galaxies is gleaned through comparing the 
variation in the total radio/UV flux ratio measurements 
of each galaxy with that derived locally for the regions of 
intense far-IR emission pinpointed by the MERLIN/VLA 
morphology. The monochromatic rest-frame 2000A lumi- 
nosities (Luv) of our SMGs are estimated from linear in- 
terpolation between the B- and i?-band magnitudes (Ta- 
ble 1). 

We estimate Lfir using the measured radio flux (Ta- 
ble 1), K-corrected using synchrotron slope of a = —0.75 
{Sy oc z/^) to restframe 1.4 GHz at the observed redshift, 
and then transformed to the far-IR using the local far-IR- 
radio correlation from Helou et al. (1985). We take this 
route to estimate the far-IR emission on arcsecond-scales 
as our submm data lack both the spatial resolution and 
full spectral information needed to estimate Lfir more di- 
rectly. We note that Garrett (2002) and Kovacs et al. 
(2004) have shown that the far-IR-radion relation does 
not seem to be change substantially out to 2: ~ 1-2.5, and 
so the approach we have adopted should be reliable. In 
this manner we predict a median total, far-IR luminosity 
of (3.5 ± 1.0) X lO^^L© for the 12 SMGs, equivalent to star 
formation rates of ^ ITOOMQyr"^ for stars more massive 
than 0.1 M0 based on a Salpeter IMF. 
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To investigate the variation of far-IR/UV ratios within 
the SMGs we first calculate the ratios for the total emis- 
sion from the galaxies. We then measure the same ratio in 
a region encompassing the peak of the radio emission us- 
ing a fixed circular aperature enclosing all radio emission 
above the 3<j contour shown in Fig. 1. We compare the to- 
tal and locally derived Lpm/Luv ratios for the 12 sources 
in Fig. 3. We see that when the regions pinpointed as the 
sites of strong activity by the MERLIN/ VL A radio mor- 
phology are considered, the obscuration levels increase by 
factors of 2-8 over the ratios derived for the whole galaxy. 
Even measured over the entire extent of the systems, the 
obscurations we derive are considerably higher than those 
seen in high redshift, restframe UV-selected populations 
which are typically Lpm/Luv ~0. 1-100 (e.g. the Lyman- 
break galaxies, Adelberger & Steidel 2000; or the BX/BM 
galaxies, Steidel et al. 2004). 
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Fig. 3. — The ratio of far-IR luminosity to restframe UV lumi- 
nosity (2000A monochromatic luminosity) plotted against total 
luminosity (far-IR + UV) for the 12 SMGs which are detected 
in the MERLIN+VLA map. We show data points correspond- 
ing to the ratio derived from total flux measurements of each 
source using MERLIN+VLA for radio and HST for the rest- 
frame UV, and based on small aperture measurements in the 
UV and radio based on the positions of the peak of the radio 
emission in the MERLIN+VLA map (Muxlow et al. 2004). 



3.1. Searching for very extended radio emission 

Having discovered that the radio emission in over half 
the SMGs in our sample is extended on V scales we now 
wish to test whether this emission extends to even larger 
scales, > 5'' or ~ 50kpc. We can constrain the fraction of 
radio flux arising on > 1.5 scales by a simple comparison 
of the 1.4-GHz fluxes of SMGs in maps from the VLA in its 
B-array (5'' synthesized beam) and A- array (1.5'' beam) 
configurations. The higher resolution A- array observations 



may resolve out a portion of the total emitted flux den- 
sity of any source which has significant extended structure 
on scales larger than the 1.5'' VLA A- array beam. As no 
B-array coverage is available for the HDF-N we have as- 
sembled VLA B-array and A- array fluxes for a sample of 
50 /iJy radio sources, including 5 submm galaxies, from 
the our SSA22 survey field (see Chapman et al. 2004). 

We calculate the median and bootstrap errors on the 
flux ratios for the radio sources from the data taken at the 
two resolutions. We start by confirming that sources show 
either similar or higher fluxes in the B-array map com- 
pared to the A-array data, as expected considering the 
median radio angular sizes of /iJy sources (Windhorst et 
al. 1993) We find median ratios of <6'b/5'a>= 1.45 ±0.27 
and that for 30% of the sample agree within 2-(j, suggest- 
ing there is no calibration offset between the fluxes from 
the two maps. For the five submm-detected sources we 
obtain a median flux ratio of < Sb/Sa >= 1-22 ± 0.12, 
compared to < Sb/Sa >= 1-22 ± 0.11 for an i?-band 
matched sample of 8 radio sources undetected or unob- 
served in the submm waveband. This suggests that the 
submm-luminous section of the /iJy radio populations is 
no more extended than the general population, after we 
have removed low redshift galaxies (with bright apparent 
magnitudes). Moreover, we see that there is only weak 
evidence that the radio emission in the SMGs is extended 
on ^ lOkpc scales. 

4. DISCUSSION & CONCLUSIONS 

The far-IR emission from local ULIRGs generally arises 
from the nuclear regions of the systems, while the UV 
light is more extended, although it contributes a negligi- 
ble fraction of the bolometric emission (Goldader et al. 
2002; Surace & Sanders 2000). We have been able to trace 
the distribution of these two components in similarly lu- 
minous galaxies at high redshifts through the combination 
of sub-arcsecond imagery in the radio and restframe UV 
from our MERLIN/ VLA and HST observations. These 
maps allowed us to investigate the relative distribution 
of obscured and unobscured star formation on kpc-scales 
within luminous submm galaxies at z ~ 2.2, when the 
Universe was only a fifth of its current age. 

As expected the UV morphologies of this sample are 
similar to the (overlapping) sample of submm galaxies im- 
aged with iJ^T/STIS and analysed by Chapman et al. 
(2003b), as weh as the ACS imaging in Small et al. (2004). 
The sample exhibits irregular and frequently highly com- 
plex morphologies in their restframe 2000A emission 
compared to optically-selected galaxies at similar redshifts, 
and have scale lengths far in excess of comparably lumi- 
nous local galaxies. 

Turning to the radio morphologies, we find that in 
- 70% (8/12) the MERLIN/VLA radio exhibits resolved 
radio emission on ~0.5-l" scales (~ lOkpc) which mirrors 
the general form of the restframe UV morphology seen 
by HST. We interpret this as strong support for the ra- 
dio emission tracing spatially-extended, massive star for- 
mation within these galaxies. This situation is very un- 
like local ULIRGs, where the high surface brightness far- 
IR/radio emission is restricted to a compact nuclear re- 
gion with an extent of less than ~ 1 kpc (Charmandaris et 
al. 2002). A more detailed analysis of the distribution of 
UV and radio emission within these galaxies (Fig. 2 & 3) 
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shows that the correspondence is rarely one-to-one, with 
variations of the UV/radio flux ratio of factors of a few on 
kpc-scales within galaxies. It is likely that comparisons us- 
ing even higher resolution data would show even stronger 
variations, as are seen in local ULIRGs (e.g. Bushouse et 
al. 2002), but which are diluted at the current resolution. 
Similarly, there is only weak evidence for correlations be- 
tween the restframe UV colors and the positions of the 
radio emitting regions within these galaxies. Here, longer 
wavelength near-IR observations with iJS'T/NICMOS, or 
from the ground, may reveal variations in the UV spectral 
slope that would show a better correlation of reddening 
with radio intensity (Small et al. 2004). 

In the remaining ^30% (4/12) of SMGs, the radio emis- 
sion is much more compact and is essentially unresolved 
- suggesting it arises in a region with a scale size of or- 
der ~ Ikpc or less (Fig. 1). In two of these cases, the 
compact radio emission is centered on a bright UV source 
(in one case clearly the nucleus of a face-on spiral galaxy, 
which is a strong X-ray source and is also the only one 
of the sources in our sample which shows AGN signatures 
in its UV spectrum), while in the other two systems the 
compact radio component is spatially offset by several kpc 
from the UV source. These configurations reflect either 
compact, nuclear starbursts and/or a dominant contribu- 
tion from an AGN to the radio emission. In half of these 
cases the AGN/nuclear starburst is also strongly obscured 
at restframe wavelengths of ~ 2000A. 

Stevens et al. (2003, 2004) have recently presented ev- 
idence for submm emission resolved on ~ 100 kpc scales 
(including apparent filaments) in the rare and extreme en- 
vironments around powerful radio galaxies and absorbed 
QSOs at z ~2-4. Our results demonstrate that in many 
cases (^70%) the far-IR emission (as seen in our MER- 
LIN/VLA radio maps) of the general SMG population is 
extended on scales ~ 10 kpc. However, the interferometric 
measurements are not suited to measuring larger scale, dif- 
fuse emission. Our VLA B-array versus A-array compar- 
ison (§3.2), however, suggests that the typical field SMG 
does not have submm emission (as traced by the radio) 
extending on scales much larger than >1''(> 10 kpc). 

Our high resolution radio and optical imaging allows 
us to address the relative obscuration of the galaxies. 
Adelberger & Steidel (2000) have suggested that high- 
luminosity galaxies at high redshift have much stonger 
obscuration than lower-luminosity galaxies, as measured 
by the ratio of far-IR to restframe-UV luminosity. This is 
certainly true on large-scales in the submm galaxies. On 
smaller scales within the SMGs, we have seen that the ob- 
scuration is roughly ~ 2 x higher over the region of intense 
radio (and by implication far-IR) emission, compared to 
the average over the whole galaxy. This suggests that there 
is highly structured reddening within the submm galaxies, 
such anisotropic obscuration would be a natural conse- 
quence of channels being blown through the dust around 
the star-formation regions by vigorous winds. 

It is worth considering that few star-burst galaxies with 
Lfir 4 X 10^^ L0 exist locally; most galaxies in our 
neighborhood with these luminosities have strong and ob- 
vious AGN components. However, at z ^ 2.2, the median 



redshift of the radio-selected SMGs, the most active galax- 
ies were evidently forming stars at rates of ^ITOOM© yr~^ 
in regions extending over ^40kpc^. The large physical 
extent of this activity contrasts markedly with the com- 
pact, nuclear starbursts typical of local redshift ULIRGs. 
This suggests that the some of the observational prop- 
erties of the star formation activity in these galaxies 
(e.g. mix of dust temperatures, ease of superwind gen- 
eration, etc.) may differ markedly from that seen in local 
"analogs". However, we also note that the star forma- 
tion surface density inferred from our radio observations 
is ~45 M0 yr~^ kpc~^, comparable to the upper-limit es- 
timated for such activity in local starburst galaxies by 
Meurer et al. (1997). This argues that the small-scale 
physical mechanisms which limit the star formation pro- 
cess within these galaxies are similar to those operating in 
the most vigorous systems locally. 

While evidence for massive amounts of molecular gas in 
submm galaxies has now been established (Frayer et al. 
1998; Neri et al. 2003; Greve et al. in preparation), and 
X-ray luminosities are consistent with a dominant role for 
star formation in the energetics of SMGs (Alexander et al. 
2004), our discovery of spatially extended radio morpholo- 
gies is perhaps the strongest piece of evidence that star 
formation dominates the bolometric output of the major- 
ity of the submm galaxy population. 

In summary, we have compared the restframe UV and 
radio morphologies on sub-arcsecond scales of a small sam- 
ple of highly luminous, dusty galaxies for which precise 
redshift s are available. This analysis shows that the radio 
emission, which we adopt as a proxy for the far-IR emis- 
sion, in these galaxies is resolved in the majority of galaxies 
- implying that dust heating (and by implication, massive 
star formation) is occuring on ~ 10 kpc scales within these 
systems. Currently, this represents our only constraint on 
the likely submm morphology of these galaxies, and one 
which will not be further testable until ALMA comes on- 
line. The overall structure of the radio emission matches 
that seen in the restframe UV - although there are strong 
variations in the relative emission on kpc-scales - which 
we interpret as resulting from highly structured dust ob- 
scuration within the galaxies. This structured obscuration 
may reflect from anisotropic dispersal of the dust as super- 
winds driven by the star formation activity blows channels 
through the inter galactic medium. Such channels would 
provide the opportunity for Lyce photons to escape from 
these otherwise highly-obscured systems, explaining the 
unexpected strength of this line in their spectra (Chap- 
man et al. 2003a, 2004). 
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Table 1 

Properties of the submm galaxies 

" The radio source extends over 3.3^^ predominantly to the west of the galaxy. Some very low surface brightness emission is 
missed in the MERLIN image which just shows three components. There is extended emission between these sources. 
^ Two resolved radio sources trace the two faint UV-detected sources within a. 1" region. 

^ Two radio sources are found within the SCUBA error circle, the brighter one listed in the table lacks an optical ID, 

and the fainter one associated with a linear galaxy 2" to the east (see Fig. 1). 

^ This source was identified in Hughes et al. (1998) as HDF 850.2 

^ HST imaging from STIS, otherwise with ACS. 

f Radio morphologies: extended (E) or compact (C). 



Source 


z 


B 


R 


Ssso 


Si. 4 


Morphology-^ 










(mJy) 


(MJy) 


Radio 


Optical 


SMM J123606.9+621021 


2.509 


25.6 


25.2 


11.6zb3.5 


74.4zb4.1 


E 


Disturbed, merger 


SMM J123616.2+621514^ 


2.578 


26.8 


25.7 


5.8zbl.l 


53.9zb8.4 


E 


3 components 


SMM J123622.7+621630^ 


2.466 


25.6 


25.4 


7.7±1.3 


70.9±8.7 


E 


Merging disks 


SMM J123629. 1+621046'' 


1.013 


26.1 


24.6 


5.0±1.3 


81.4±8.7 


E 


Disturbed 


SMM J123655.8+621200^ 


2.743 


25.4 


25.3 


8.0±1.8 


21.0±6.2 


E 


Disturbed, merger 


SMM J123707.2+621408 


2.484 


26.9 


26.0 


6.3zbl.3 


45.3zb7.9 


E 


Blue & red pair 


SMMJ123712.0+621325 


1.992 


26.0 


25.8 


4.1zbl.3 


53.9zb8.1 


E 


Disturbed with dusty component 


SMM J123712.1+621212^ 


2.914 


27.0 


25.5 


8.0±1.8 


21.0±4.0 


E 


Disturbed, double source 


SMM J123618. 3+621551^ 


1.865 


26.0 


25.9 


7.3±1.1 


150.5±11.2 


C 


Small group 


SMM J123621.3+621708^'^ 


1.988 


25.1 


24.9 


7.8±1.9 


148.1±11.2 


C 


Linear 


SMM J123635.6+621424 


2.005 


24.2 


24.2 


5.5zbl.4 


87.8zb8.8 


C 


Disturbed, face-on spiral 


SMM J123646.1+621449 


1.7 


25.8 


25.7 


10.3=b2.2 


124.3zb7.9 


C 


Disturbed (photo-z) 
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